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Abstract 

In the context of the minimal supersymmetric extension of the Standard Model 
the effect of a realistic wall profile is studied. It has been recently showed that in 
the presence of light stops the electroweak scale phase transition can be strong 
enough for baryogenesis. In the presence of non-trivial CP-violating phases of 
left-handed mixing terms and Higgsino mass, the largest ub/s is created when 
Higgsino and gaugino mass parameters are degenerate, [i = M^. In the present 
paper we show that realistic wall profiles suppress the generated baryon number 
of the universe, so that quite a stringent bound | sin (j)^ > 0.2 for /x-phase 
can be inferred. 
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The Minimal Supersymmetric Standard Model (MSSM) has appeared to be one 
of the most promising candidates to explain the observed baryon asymmetry of the 
universe ub/s ~ 10~ 10 generated at electroweak scale. [l|]. Although all requirements 
are included already in the Standard Model M, [J, the phase transition 0] has appeared 
to be too weakly first order to preserve the generated baryon asymmetry ||. Also it has 
been shown that the CP- violation needed for baryogenesis is too small in the Standard 
Model pi. Therefore some new physics besides the Standard Model is necessarily 
needed, provided that the baryon asymmetry is generated during the electroweak 
phase transition. 

Because MSSM is one of the most appealing extensions of the Standard Model, 
it has been worthwhile to study whether it is possible to generate and preserve the 
baryon asymmetry in it. Indeed, recent analyses show that there exists a region of 
the parameter space where the phase transition is strong enough |3j. It is required 
that tan/5 < 3, the lightest stop is lighter than the top quark and the lightest Higgs 
must be detectable by LEP2. The bound given above may relax due to two and higher 
-loop effects, which seem to strengthen the phase transition H. Unlike the Standard 
Model where the source of CP-violation is solely the Cabibbo-Kobayashi-Maskawa 
matrix, MSSM contains an additional source due to the soft supersymmetry breaking 
parameters which are related to the stop mixing angle. 

In a recent paper Carena et al. || has analysed the region of supersymmetric 
parameter space where the baryon asymmetry generated at the electroweak scale is 
consistent with the observed one. The generation and survival of large enough %/s ~ 
4 x 1CT 11 seem to require that Mz < mi < m t , the mass of the lightest Higgs boson 
is bounded by tuh < 80 GeV whereas CP-odd boson has mass m A > 150 GeV. The 
analysis gave dependence of hb/s on Higgsino mass parameters and its phase M 
as well as gaugino mass parameters Mi and M 2 . The optimal choice of parameters 
showed up to be \fi\ ~ M 2 , | sin0 M | > 0.06, and the dependence on Mi is weak, so 
that it can be chosen to be equal to M 2 . In the analysis of Carena et al. it was 
chosen the left-handed stop mass parameter rriQ to be mq = 500 GeV, effective soft 
supersymmetry breaking parameter A t = and right-handed stop mass parameter 
mu = —fhjj < 



to obtain the most optimistic bounds (i.e. maximize Ub/s) [J With these parame- 
ter values the CP-violating source is generated essentially by Higgsino and gaugino 
currents and the right-handed stop contribution is negligible, so that without a loss 
of generality one can set sin(0 M + <p A ) = 0. The bound ([!]) is due to the colour 
non-breaking condition, i.e no colour breaking minimum must be deeper than the 

1 About restrictions and validity of these results, see ||. 
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normal electroweak breaking (and colour conserving) minimum. It is defined at zero- 
temperature, thus vq = 246.22 GeV. Also it was used v w = 0.1 and L w = 25/T. With 
these conditions a large enough baryon asymmetry could have been created. 

In the paper || the baryon asymmetry was inferred by first calculating the CP- 
violating sources and then solving the relevant Boltzmann equations. It was shown 
that the baryon to entropy ratio reads 

— = -g{h) — 2 — ' ( 2 ) 

s vis 

where g(ki) is a numerical coefficient depending on the degrees of freedom, D the 
effective diffusion rate, T ws = Skol^T (k = 1) || the weak sphaleron rate, v w the wall 
velocity. The entropy density s is given by 

2vr 2 ^T 3 

° = —45—. (3) 

where g* s is the effective number of relativistic degrees of freedom. The coefficient A 
was shown to be a certain integral over the source j(u) = v w f(ki)d u J°(u): 

A= ~mT duj(u)e- x + u , (4) 
D\ + Jo 



where A + = (v w + yv^ + ATD)/(2D) and the wall was defined so that it begins at 
u = 0. Here u is the co-moving coordinate u = z+v w t supposing that the wall moves in 
the direction of z-axis. (f(ki) is again a coefficient depending on the number of degrees 
of freedom present in thermal path and related to the definition of the effective source 
[§, 0.) Thus the coefficient A is dependent on the actual wall shape via 



where H = \J Hf + iff? tan/3 = H1/H2 and ifj's are the real parts of the neutral 
components of the Higgs doublets. In || the wall shape was, however, taken ad hoc. 
It was assumed to have a simple sinusoidal form so that the field H(u) can be given 
I'.v 

[9(u)-e{u-L w )}+v9(u-L w ) (6) 



H s (u) = \ 



1 — cos 



un 



and the field angle f3{u) by 



AM = ¥ 



1 — COS 



U71 



[6{u)-8{u-L w )] + Al36{u-L w ), (7) 



where A/3 is given by A/? = (3 (To) — arctan(mi(To)/m 2 (T )), calculated at the tem- 
perature where curvature of the one-loop effective potential vanishes at the origin. 



Inserting these Anzdtse to Eq. ([5]), we obtain a corresponding contribution to the 
CP-violation 

Using this approximation, the results of || was inferred. 

In the present paper a more realistic prescription of wall shape is used. Working at 
the critical temperature T c and using the one-loop resummed effective potential [|ll|] we 
numerically find the path of smallest gradient 7 5 from (Hi, H 2 ) = (0, 0) to (Hi, H 2 ) = 
(vi,v 2 ) = (v i(T c ), t> 2 (T c )), which well approximates the true solution. Moreover, using 
path 7 S , a upper bound for / is necessarily obtained. The true solution lies necessarily 
between 7 9 and straight line from (0,0) to (vi,v 2 ) (which leads to n B = 0) as can be 
concluded by studying the LagrangeanQ. Thus |/| along such a path is smaller than 
along j g . 

The effective potential for MSSM at finite temperature can be expressed in three 
parts 

V eff (<j>, T) = V (<P) + + V ltT {<f>) + AVt{4>), (9) 

where Vq is the tree level zero-temperature potential, V\ the renormalized 1-loop zero- 
temperature potential, V\ 7 t the 1-loop finite temperature potential and AVr the daisy- 
resummed part. They are given by 

2 /2 

V (H) = m\H\ + m\H\ + 2m\ 2 H x H 2 + 9 —^-(H\ - H 2 2 f, (10) 

o 

t,ii,2,W,Z 



and 



A^t(H) = -^I>K 3 (lr,T) -m?(#)], (13) 



where m^H) and rrii(H,T) are the zero temperature and temperature corrected field 
dependent masses, respectively, rij are the degrees of freedom of each particle (including 
— sign for fermions) and 



J.( x 2 ) = / dyy 2 \n(l±e-^+* 5 ). (14) 



Note, that we have neglected the b-quark, b-squark as wall as other generation contri- 
butions as small ones. The heavy supersymmetric particles do not contribute neither. 

2 This is true providing that the field is not strongly oscillating within the wall, but is a smooth 
configuration. 



The mass parameters mi, m.2, m.12, are related to /3, rriA, itih and other parameters 



of the theory, as given in ||11|| . Using these formulas, we can solve 7 g and define the 



corresponding CP-violation integral 

A = f du L {H{u)2 ^r )e ~ x+u (15) 

at the critical temperature along the path j g . It also shows up that the form of 
profiles is, with good acccuracy, the form of a kink. Indeed, if we reparametrise the field 
(Hi, H 2 ) to a component pointing towards (vi, t> 2 ), Hn and to a component orthogonal 



to that, H±, it appears that the ratio of maximum value of H± to v — \lv\ + v\ is in 
any case smaller than 0.004. Thus the bending of the path, I.e. the deviation from a 
straight line is small. (Note, that here v is the value of vacuum at T c , thus not equal 
to v = 247 GeV.) 

The ratio 7 7 / I s gives immediately the supression of n B /s with respect to the results 
of Carena et al. ||. Hence the value needed for soft supersymmetry mixing phase 
sin^ is increased by factor I 7 /I a . In the Figure 1. we have presented the value of 
the integral I s as a function of fi for several values of tua- From the figure it can be 
read out that increasing /x decreases the (absolute) value of I 7 so that for /j > 250 
GeV inequality |I 7 | < 1 holds. Also it can be found that increasing rriA with factor r 
decreases J 7 roughly by factor 1/r within the parameter range studied. This behaviour 
is likely be more general than just restricted to analysis of ||, because the amount of 
CP-violation is in general proportional to the change of j3 on the path from the origin 
to the non-trivial vacuum. 

In Figure 2. a comparison to the result of Carena et al. || is made. We have 
plotted the ratio I 7 /I s also as a function of \x with several values of tua- There appears 
a clear tendency of additional suppression: for very small /1, I-y/I s ~ 0.4, whereas for 
/j ~ 250 GeV, I y /I a ~ 0.1. For optimal values of //, 150 GeV < /i < 250 GeV the 
supression factor is at least 0.3. Also it appers that the dependence of this suppression 
factor on CP-odd boson mass rriA is relatively weak. Taking in to the account that 
according to the analysis of || it is required that | sin0^| > 0.06, this is now converted 
to a more stringent bound | sin^l > 0.2 which remarkably weakens the possibility of 
baryogenesis in MSSM. However, for fi ~ 250 GeV the bound rises to | sin0^| > 0.6 
which possibly is already too large. 

In the present paper we have calculated the amount of CP-violation in the bubble 
wall of the minimal supersymmetric extension of the Standard Model at the elec- 
troweak scale phase transition. It has shown up that the previous estimates tend to 
be too optimistic and an extra suppression of (at least) 0.3 is found. This tends to 
make the electroweak baryogenesis in MSSM more difficult and less likely. However, 
more analysis on the model is needed, in particular to clarify how higher corrections 
(more scalar insertions) to the CP-violating source behave. If the higher corrections to 



the source expand in the powers of (H\dH2 — H2dH\) no help from them is expected. 
If they, however, order by some other expansion parameter, their contribution to the 
source may be remarkably large. Unfortunately this may also lead to the situation 
where non-perturbative effects are important. Also the effect of higher order corrected 
effective potential remains to be studied. The two-loop corrections tend to strenghten 
the phase transition and thus relax the bounds The two-loop contributions may, 
however, affect also directly the value of the integral / needed in the calculation of 
CP-violating source. On the other hand, also the sphaleron rate r„, s is under discus- 
sion [0 and changes on that may change significantly the conclusions made about 
baryogenesis. 
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Figure captions 



Figure 1. Values of the integral L f = j^(H 2 f3')'e x+u (where comma stands for u 
-derivative) as a function of the soft supersymmetry parameter \i = Mi = M 2 with 
v w = 0.1, m Q = 500 GeV, rh v = m^ rit , tan (3 = 2 and A t = 0. 

Figure 2. Values of the ratio I-,/I s as a function of fi = Mi = M 2 with v w = 0.1, 
m Q = 500 GeV, = m^ rit , tan/3 = 2 and A t = 0. For I s wall width L w = 25/T is 
used. 
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